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B.  TECHNICAL  SECTION 

I.  Abstract 

The  overall  objective  of  this  proposal  is  to  investi¬ 
gate  the  robustness  to  departures  from  independence  of  methods 
currently  in  use  in  reliability  studies  when  competing  failure 
modes  or  competing  causes  of  failure  associated  with  a  single 
mode  are  present  in  a  series  system.  The  first  specific  aim 
is  to  examine  the  error  one  makes  in  modeling  a  series  system 
by  a  model  which  assumes  statistically  independent  component 
lifetimes  when  in  fact  the  component  lifetimes  follow  some 
multivariate  distribution.  The  second  specific  aim  is  to  assess 
the  effects  of  the  independence  assumption  on  the  error  in 
estimating  component  parameters  from  life  tests  on  series 
systems.  In  both  cases,  estimates  of  such  errors  will  be  deter¬ 
mined  via  mathematical  analysis  and  computer  simulations  for 
several  prominent  multivariate  distributions.  A  graphical  dis¬ 
play  of  the  errors  for  representative  distributions  will  be 
made  available  to  researchers  who  wish  to  assess  the  possible 
erroneous  assumption  of  independent  competing  risks.  A  third 
aim  is  to  tighten  the  bounds  on  estimates  of  component  relia¬ 
bility  when  the  risks  belong  to  a  general  dependence  class  of 
distributions  (for  example,  positive  quadrant  dependence,  posi¬ 
tive  regression  dependence,  etc.) .  Major  decisions  involving 
reliability  studies,  based  on  competing  risk  methodology,  have 
been  made  in  the  past  and  will  continue  to  be  made  in  the  future. 
This  study  will  provide  the  user  of  such  techniques  with  a 
clearer  understanding  of  the  robustness  of  the  analyses  to  de¬ 
partures  from  independent  risks,  an  assumption  commonly  made 
by  the  methods  currently  in  use. 
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II.  Specific  Objectives 

The  overall  objective  is  to  investigate  the  robust¬ 
ness  to  departures  from  independence  to  methods  currently  in 
use  in  reliability  studies  when  competing  failure  modes  or  com¬ 
peting  causes  of  failure  associated  with  a  single  mode  are 
present  in  a  series  system.  We  shall  also  refer  to  such  com¬ 
petitive  events  as  competing  risks.  The  approach  will  be  through 
the  investigation  of  certain  aspects  of  specific  parametric  multi¬ 
variate  distributions  or  by  classes  of  distributions  which  are 
appropriate  in  reliability  analyses  when  there  are  competing 
risks  present. 

The  specific  objectives  are: 

1)  to  assess  the  error  incurred  in  modeling  system 
life  in  a  series  system  assumed  to  have  indepen¬ 
dent  component  lifetimes  when  in  fact  the  com¬ 
ponent  lifetimes  are  dependent. 

2)  to  assess  the  error  in  estimating  component  param¬ 
eters  (i.e.,  component  reliability,  mean  com¬ 
ponent  life,  etc.)  in  a  series  system  employing 
either  parametric  or  nonparametric  models  which 
assume  independent  component  failure  times  when 

in  fact  the  lifetimes  are  dependent  and  follow 
some  plausible  multivariate  distribution.* 

3)  to  derive  bounds  on  component  reliability  when 
the  failure  modes  are  dependent  and  fall  in  a 
particular  dependence  class  (e.g.,  positive  quad¬ 
rant  dependence,  positive  regression  dependence, 
etc. ) . 

4)  to  develop  tests  of  independence,  based  on  data 
collected  from  series  systems,  by  making  some 
restrictive  assumption  about  the  structure  of  the 
systems . ** 

*  A  plausible  parametric  multivariate  distribution  will  be 
one  that  satisfies  one  of  the  following  conditions: 

i)  the  distribution  of  the  minimum  of  the  component 
failure  times  closely  approximates  widely  accept¬ 
ed  families  of  system  life  distributions. 

or  ii)  the  marginal  distributions  closely  approximate 
the  distributions  of  component  failure  times  in 
the  absence  of  other  failure  modes. 

**This  objective  has  been  added  to  the  original  objectives  be¬ 
cause  it  answers  a  natural  question  raised  by  our  preliminary 
investigation. 
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III.  Introduction  to  Problem  and  Significance  of  Study 

Alvin  Weinberg  (1978)  in  an  editorial  comment  in  the 
published  proceedings  of  a  workshop  on  Environmental  Biologi¬ 
cal  Hazards  and  Competing  Risks  noted  that  "the  question  of 
competing  risks  will  not  quietly  go  away:  corrections  for  com¬ 
peting  risks  should  be  applied  routinely  to  data."  The  problem 
of  competing  risks  commonly  arises  in  a  wide  range  of  experi¬ 
mental  situations.  Although  we  shall  confine  our  attention 
in  the  following  discussion  to  those  situations  involving 
series  systems  in  which  competing  failure  modes  or  competing 
causes  of  failure  associated  with  a  single  mode  are  present, 
it  is  certainly  true  that  we  might  just  as  easily  speak  of 
clinical  trials,  animal  experiments,  or  other  medical  and  bio¬ 
logical  studies  where  competing  events  interrupt  our  study  of 
the  main  event  of  interest  (cf .  Lagokos  (1979))  . 

Consider  electronic  or  mechanical  systems,  such  as 
satellite  transmission  equipment,  computers,  aircraft,  missiles 
and  other  weaponry  consisting  of  several  components  in  series. 
Usually  each  component  will  have  a  random  life  length  and  the 
life  of  the  entire  system  will  end  with  the  failure  of  the 
shortest  lived  component.  We  will  examine  two  situations  more 
closely  in  which  competing  risks  play  a  vital  role. 

First,  suppose  we  are  attempting  to  evaluate  system  life 
from  knowledge  of  the  individual  component  lifetimes.  Such 
an  evaluation  will  utilize  either  am  analysis  involving  math¬ 
ematical  statistics  or  a  computer  simulation.  At  a  recent 
conference  on  Modeling  and  Simulation,  McLean  (1981)  presented 
a  scheme  to  simulate  the  life  of  a -missile  which  consisted  of 
many  major  components  in  series.  The  failure  distribution  asso¬ 
ciated  with  each  component  was  assumed  to  be  known  (usually 
exponential  or  Weibull'. )  To  arrive  at  the  system  failure  dis¬ 
tribution,  the  components  were  assumed  to  act  independently  of 
each  other.  Realistically,  this  may  or  may  not  be  the  case. 

If  the  component  lifetimes  were  dependent  for  any  reason,  the 
computed  system  failure  distribution  (as  well  as  its  subsequent 
parameters  such  as  system  mean  life  and  system  reliability  for 
a  specified  time)  would  only  crudely  approximate  the  true 
distribution.  The  first  specific  aim  of  this  proposal  is  to 
ascertain  the  error  incurred  in  modeling  system  life  in  a 
series  system  assumed  to  have  independent  component  lifetimes 
(i.e.,  risks)  when,  in  fact,  the  risks  are  dependent. 

Second,  suppose  we  wish  to  evaluate  some  aspect  of  the 
distribution  of  a  particular  failure  mode  based  on  a  typical 
life  test  of  a  series  system.  The  response  of  interest  is  the 
time  until  failure  of  a  particular  mode  of  interest.  Frequently 
this  response  will  hot  be  observable  due  to  the  occurrence  of 
some  other  event  which  precludes  failure  associated  with  the 
node  of  interest.  We  shall  term  such  competing  events  which 
interrupt  our  study  of  the  main  failure  modes  of  interest  as 
competing  risks . 
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Competing  risks  arise  in  such  reliability  studies  when 

1)  the  study  is  terminated  due  to  a  lack  of  funds  or  the 
pre-determined  period  of  observation  has  expired 
(Type  I  censoring) . 

2)  the  study  is  terminated  due  to  a  pre-determined  number 
of  failures  of  the  particular  failure  mode  of  interest 
being  observed  (Type  II  censoring) . 

3)  some  systems  fail  because  components  other  than  the 
one  of  interest  malfunction. 

4)  the  component  of  interest  fails  from  some  cause  other 
than  the  one  of  interest. 

In  all  four  situations,  one  may  think  of  the  main  event  of 
interest  as  being  censored,  i.e.,  not  fully  observable.  In  the 
first  two  situations ,  the  time  to  occurrence  of  the  event  of 
interest  should  be  independent  of  the  censoring  mechanism.  In 
such  instances,  the  methodology  for  estimating  relevant  reliabili¬ 
ty  probabilities  has  received  considerable  attention  (cf.  David 
and  Moeschberger  (1978) ,  Kalbfeish  and  Prentice  (1980) ,  Elandt- 
Johnson  and  Johnson  (1980) ,  Mann,  Schafer,  Singpurwalla  (1974) 
and  Barlow  and  Proschan  (1975)  for  references  and  discussion)* 

In  the  third  situation,  the  time  to  failure  of  the  component  of 
interest  may  or  may  not  be  independent  of  the  failure  times  of 
-other  components  in  the  system.  For  example ,  there  may  be 
common  environmental  factors  such  as  extreme  temperature  which 
may  affect  the  lifetime  of  several  components.  Thus  the  question 
of  dependent  competing  risks  is  raised.  A  similar  observation 
may  be  made  with  respect  to  the.  fourth  situation,  viz.,  failure 
times  associated  with  different  failure  modes  of  a  single  com¬ 
ponent  may  be  dependent.  For  a  very  special  type  of  dependence, 
the  models  discussed  by  Marshall-Olkin  (1967) ,  Langberg,  Proschan 
and  Quinzy  (1978) ,  and  Langberg,  Proschan,  and  Quinzy  (1981) 
allow  one  to  convert  dependent  models  into  independent  ones. 

If  no  assumptions  whatever  are  made  about  the  type  of  * 
dependence  between  the  distribution  of  potential  failure  times, 
there  appears  to  be  little  hope  of  estimating  relevant  component 
parameters.  In  some  situations,  one  may  be  appreciably  misled 
(cf.  Tsiatis  (1975),  Peterson  (1976)).  However,  as  Easterling 
(1980)  so  clearly  points  out  in  his  review  of  Bimbaum's  (1979) 
stenograph 

"there  seems  to  be  a  need  for  some  robustness 
studies.  How  far  might  one  be  off,  quantita¬ 
tively,  if  his  analysis  is  based  on  incorrect 
assumptions?” 

The  second  specific  aim  will  address  this  important 
issue.  First  if  a  specific  parametric  model  which  assumes 
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independent  risks  has  been  used  in  the  analysis,  it  would 
be  of  interest  to  know  how  the  error  in  estimation  is 
affected  by  this  assumption  of  independence.  That  is,  if 
independent  specific  parametric  distributions  are  assumed 
for  the  failure  times  associated  with  different  failure 
modes  when  we  really  should  use  a  bivariate  (or  multivariate) 
distribution,  then  what  is  the .magnitude  of  the  error  in 
estimating  component  parameters?  Secondly,  one  may  wish  to 
allow  for  a  less  stringent  type  of  model  assumption,  and  ask 
the  same  question  with  regard  to  the  estimation  error.  That 
is,  if  a  nanparametrlc  analysis  is  performed,  assuming  in¬ 
dependent  risks,  when  some  types  of  dependencies  may  be 
present,  then  what  is  the  magnitude  of  the  estimation  error? 

The  third  specific  aim  will  attempt  to  obtain  bounds  on 
the  component  reliability  when  the  failure  times  belong  to 
a  broad  dependence  class  (e.g.,  association,  positive  quadrant 
dependence,  positive  regression  dependence,  etc.).  More 
details  will  be  presented  in  the  methods  section. 

In  summary,  competing  risk  analyses  have  been  performed 
in  the  past  and  will  continue  to  be  performed  in  the  future. 
This  study  will  provide  the  user  of  such  techniques  with  a 
clearer  understanding  of  the  robustness  to  departures  from 
independent  risks,  an  assumption  which  most  of  the  methods 
currently  in  use  assume. 
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IV.  Progress  Report  on  First  Year's  Work 

We  believe  that  substantial  progress  has  been  made  in 
dealing  with  the  first  two  objectives  as  outlined  on  page  4.  We 
have  examined  the  consequences  of  erroneously  assuming  independence 
when  the  true  multivariate  distribution  of  component  lifetimes 
follows  either  a  Gumbel  (1960)  bivariate  exponential  model  or  the 
well-known  Marshall-Olkin  (1967)  bivariate  exponential  model. 

These  bivariate  distributions  model  two  series  systems  which  have 
quite  different  failure  mechanisms.  The  former  has  the  flexibility 
of  allowing  either  positive  or  negative  correlation  (not  due  to 
simultaneous  failure)  between  two  component  lifetimes  in  a  series 
system.  The  latter  allows  dependency  to  enter  via  simultaneous 
failure  of  the  components.  Results  associated  with  both  models 
have  been  submitted  for  publication.  (See  Appendices  A  and  B  for  copies  of 
these  papers.)  The  paper  dealing  with  the  Gumbel  model  has  recently 
been  substantially  revised  for  Technometrics .  In  addition  to  clarifying 
certain  aspects  of  the  paper  as  originally  submitted,  we  have  investi¬ 
gated  how  well  the  Mann-Grubbs  (1974)  confidence  intervals  perform 
for  varying  degrees  of  correlation  among  the  components  and  for  various 
sample  sizes.  The  Gumbel  results  were  presented  in  an  invited  talk 
at  the  Applied  Statistics  Conference  in  Newark  in  December,  1982. 

The  Marshall-Olkin  results  were  presented  at  the  joint  ASA-IMS-Biometric 
Society  meetings  in  Toronto  in  August,  1983.  Also,  some  aspects  of 
these  results  were  given  in  an  invited  session  on  Survival  Analysis 
at  the  AAAS  in  Detroit  in  May,  1983. 

We  also  have  examined,  in  a  related  issue,  the  asymptotic 
bias  of  the  Kaplan-Meier  product  limit  estimator  under  dependent 
competing  risks  and  have  evaluated  the  magnitude  of  this  for  some 
specific  bivariate  distributions.  These  results  have  also  been 
submitted  for  publication.  (See  Appendices  for  a  copy.) 

Work  is  near  completion  with  respect  to  evaluating,  the 
consequences  of  erroneously  assuming  independence  when  modeling  system 
reliability  from  ccnx>lete  component,  information.  This  work  includes,  in  addition 
to  the  above  models",  Downton's  bivariate  exponential,  two  additional 
Gumbel  models,  and  the  bivariate  model  proposed  by  Oakes  (1982). 


A  natural  question,  which  arose  from  our  investigation 
into  the  effects  of  dependence  on  modeling  system  reliability  and 
in  estimating  component  parameters  from  system  data,  is  how  can 
one  test,  statistically,  the  assumption  of  independence?  A  test 
of  independence,  based  on  data  collected  from  series  systems,  is 
not  possible,  unless  some  restrictive  assumptions  are  made,  due 
to  the  nonidentifiability  problem.  That  is,  for  any  set  of  ob¬ 
servable  information  coming  from  a  series  system  with  dependent 
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component  lifetimes  there  exists  a  system  with  independent  com¬ 
ponent  lives  which  has  the  same  observable  quantities  (c.f.  Tsiatis 
(1975),  Rose  (1973),  Basu  and  Klein  (1982)).  By  making  some 
restrictive  assumptions  about  the  joint  distribution  of  the  com¬ 
ponent  lifetimes  or  about  the  net  component  lifetimes  it  may  be 
possible  to  test  for  independence  in  this  restrictive  setting. 

An  approach  which  makes  some  assumptions  about  the 
marginal  component  lifetimes  and  then  uses  the  information  to 
modify  the  standard  nonoarametric  test  for  independence  is  presented 
in  another  paper  in  the  Appendix  D.  The  thrust  of  our  future 
efforts  will  be  to  study  further  aspects  of  the  test  system 
reliabilities  under  various  conditions.  In  conclusion,  we  believe 
that  we  are  on  target  with  respect  to  dealing  with  the  objectives 
as  outlined  in  our  original  proposal  as  well  as  the  one  added 
objective  which  serves  to  enhance  the  overall  effort. 

VI .  Methods 

We  refer  to  pages  8052  of  the  original  proposal  for  a 
discussion  of  the  general  methodology. 
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Abstract 


This  paper  investigates  the  consequences  of 
departures  from  independence  when  the  component 
lifetimes  in  a  series  system  are  exponentially 
distributed.  Such  departures  are  studied  when 
the  joint  distribution  is  assumed  to  follow  a 
Gumbel  bivariate  exponential  model.  Two  distinct 
situations  are  considered.  First,  in  theoretical 
modeling  of  series  systems,  when  the  distribution 
of  the  component  lifetimes  is  assumed,  one  wishes 
to  compute  system  reliability  and  mean  system 
life.  Second,  errors  in  parametric  and  nonpara- 
metric  estimation  of  component  reliability  and 
component  mean  life  are  studied  based  on  lifetest 
data  collected  on  series  systems  when  the  assump¬ 
tion  of  independence  is  made  erroneously. 

Keywords:  competing  risks,  component  life,  modeling 

series  systems,  robustness  studies,  system 
reliability,  Gumbel  bivariate  exponential 


1.  INTRODUCTION 


Consider  a  system  consisting  of  several  components  linked 
in  series.  For  such  a  system  the  failure  of  any  one  of  the 
components  causes  the  system  to  fail.  Common  assumptions  made 
in  modeling  and  analyzing  data  from  such  a  system  are  that  the 
component  lifetimes  are  independent  and  exponentially  distri¬ 
buted.  Many  authors  have  considered  the  problem  of  analyzing 
a  series  system  with  exponential  component  lives.  For  example, 
confidence  bounds  for  system  reliability  assuming  independent 
exponentially  distributed  component  lifetimes  were  presented 
in  Mann(1974)  and  Mann  and  Grubbs  (1974).  (cf . ,  see  Mann, 

Schafer,  and  Singpurwalla  (1974)  for  a  more  comprehensive  review.) 
More  recently,  work  invoking  the  assumption  of  independent 
exponentially  distributed  lifetimes  has  been  presented  by  Chao 
(1981)  and  Miyamura  (1982) .  Estimation  of  component  parameters 
from  series  system  data  has  been  treated  by  Boardman  and  Kendell 
(1970)  in  the  context  of  independent  exponential  component  lives. 
Some  authors  suggest  a  nonparametric  alternative  to  the  estima¬ 
tion  of  component  reliability  based  on  series  system  data  (cf. 
Kalbfleisch  and  Prentice  (1980)  and  Lawless  (1982)). 

The  assumption  of  independence  is  essential  to  these  analyses 
and  also  an  important  concern.  Several  authors  have  shown  that 
this  assumption,  by  itself,  is  not  testable  since,  based  on  data 
from  a  series  system,  there  is  no  way  to  distinguish  between  an 


2. 


independent  and  dependent  model.  (See  Tsiatis  (1975),  Peterson 
(1976) ,  and  Basu  (1981)  for  a  discussion  of  nonidentifiability 
results.)  In  many  situations  one  may  be  appreciably  misled  by 
the  independence  assumption. 

Lagakos  (1979)  ,  in  a  study  of  the  effects  of  various  types 
of  dependence  among  component  lifetimes,  has  noted  that  most 
methods  of  analysis  have  assumed  noninformative  models  of  which 
independence  is  a  special  case.  He  points  out,  "it  is  important 
to  be  aware  of  the  possible  consequences  of  making  this  assumption 
when  it  is  false."  Furthermore,  Easterling  (1980)  states  in  his 
review  of  Birnbaum's  (1979)  monograph  on  competing  risks  "there 
seems  to  be  a  need  for  some  robustness  studies.  How  far  might 
one  be  off,  quantitatively,  if  his  analysis  is  based  on  incorrect 
assumptions?" 

In  this  paper  we  consider  the  consequences  of  departures  from 
independence  when  the  component  lifetimes  are  exponentially  dis¬ 
tributed.  Such  dep  irtures  from  independence  may  be  related  to 
some  common  environmental  factor,  only  present  when  the  components 
are  linked  together  in  series.  The  load  each  component  is  subject 
to  is  either  reduced  or  increased  according  to  the  age  of  the 
system.  To  study  such  departures  we  have  selected  a  model  pro¬ 
posed  by  Gumbel  (1960).  Gumbel's  model  retains  the  assumption 
of  exponentially  distributed  component  lifetimes  while  allowing 
the  flexibility  of  both  positive  or  negative  mild  correlation 
between  component  lifetimes. 
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The  effects  of  a  departure  from  the  assumption  of  independent 
component  lifetimes  in  a  series  system  will  be  addressed  in  two 
distinct  situations.  The  first  situation  arises  in  modeling  the 
performance  of  a  theoretical  series  system  constructed  from 
components  whose  lifetimes  are  exponentially  distributed.  Here, 
based  on  testing  each  component  separately  or  based  on  engineering 
design  principles,  it  is  reasonable  to  assume  that  the  components 
are  known  to  be  exponentially  distributed  with  known  parameter 
values.  Based  on  this  information  we  wish  to  calculate  parameters 
such  as  the  mean  system  life  or  system  reliability  of  a  series 
system  constructed  from  these  components.  In  Section  2  we  describe 
how  the  values  of  these  quantities  are  affected  by  departures 
from  independence  when  the  component  parameters  are  completely 
specified.  In  Section  3  we  study  the  performance  of  the  Mann- 
Grubbs  (1974)  confidence  bounds  on  system  reliability  for  small 
sample  sizes  and  for  varying  degrees  of  correlation  .when  the  component 
parameters  are  estimated  from  component  data. 

The  second  situation  involves;  making  inferences  about  compo¬ 
nent  lifetimes  distributions,  component  reliabilities,  and  compo¬ 
nent  mean  lives,  from  data  collected  on  series  systems.  Commonly, 
data  collected  on  such  systems  are  analyzed  by  assuming  a  constant- 
sum  model  of  which  independence  is  a  special  case  (cf.  Williams 
and  Lagakos  (1977)  and  Lagakos  and  Williams  (1978).  In  Section  4 
we  study  the  properties  of  the  maximum  likelihood  estimators  of 
component  parameters  calculated  under  an  assumption  of  independent 
exponential  component  lifetimes  when  the  component  lifetimes  are 
Gumbel  bivariate  exponential.  Because  of  the  usage  of  nonparametric 
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estimates  of  component  reliability,  we  also  study  the  estimation 
error  of  the  Kaplan-Meier  (1958)  estimator  when  the  assumption 
of  independence  is  made  erroneously. 


2.  MODELING  SYSTEM  RELIABILITY  FROM  COMPLETE 
COMPONENT  INFORMATION 

Consider  a  two  component  series  system  with  component  life- 
lengths  X^ ,  X2-  Suppose  that  X^  has  an  exponential  survival 
function 

Fi  (t)  =  P(Xi>t)  =  exp  (-Xit), 

X ^ ,  t>0,  i  =  1,2. 

This  assumption  is  made  on  the  basis  of  extensive  testing  of 
each  component  separately  or  ^n  knowledge  of  the  underlying 
mechanism  of  failure.  The  value  of  x^  is  assumed  known.  If 
(X^,X^)  are  independent,  then  the  time  to  system  failure  has 
an  exponential  distribution  with  failure  rate  X=X^+\2,  and 
the  system  reliability  is  given  by 

Fj  (t)  =  P  jmin  (X,  ,  X2)  >t  i  independence^  exp  (-xt).  (2.1) 

Suppose  that  the  actual  joint  distribution  of  (X^,X2)  has 
the  form  proposed  by  Gumbel  (1960)  ,  namely, 

P(X1>x1,X2>x2) 

*{exp  (-X^Xj -X2x2)*;  [i+o  (1-exp  (1-exp  (-\2x2))J  (2.2) 

and  the  joint  probability  density  of  (Xj, is 

f  (x^  ,x2)*X1X2  [exp  (-X1x1-X2x2)j  j^l+a  (2  exp  (-XjX^-l)  (2  exp  (-X2x2)-l)j  (2.3) 
where  in  both  (2.2)  and  (2.3),  x^,x2,X^,  x2>0,-l<_a<l.  This 


distribution  has  marginal  survival  functions  equivalent  to  those 
for  the  independent  model  which,  in  part,  is  the  reason  for 
choosing  it.  The  correlation  between  (X-^jX^)  is  p  =a/4  and  a=0 
is  equivalent  to  (X^,X2)  being  independent.  For  p  >0  (<0)  the 
components  are  positively  (negatively)  quadrant  dependent  (cf. 
Barlow  and  Proschan  (1975)).  Furthermore,  the  conditional 
expectation  of  X-^,  given  X2=x2,  is 

E(X1|X2=x2)=  j  Q.+2p-4pexp(-X2x2)j  . 

If  (X^,X2>  have  the  joint  distribution  (2.3)  then  the  true 
system  reliability,  is 

Fd  (t)  =  P  |jnin(X^  ,X2)  >t  (dependence  J 

=  exp (- At)  Jj.+4p  (1-exp (-x^t)  )  (l-exp(-x2t)  )”7.  (2.4) 

From  (2.1)  and  (2.4)  we  can  easily  see  the  amount  of  error 
in  modeling  system  reliability  is 
A(t)=FD(t)-FI(t) 

=4p Q-exp (-X^tjJ  Jj.-exp(-X2t)J  exp£- (  X1+x2)  tj .  (2.5) 

Note  that  |  A ( t )  |  increases  as  |p|  increases,  for  fixed  X^,X2,  and  t. 

The  magnitude  of  A(t)  ,  of  course,  depends  upon  x^rX2f  anc^  o* 

When  x1=X2='J>.  one  can  show  that  A(t)  is  maximized  at 
t  =  fin  2l/p  (fixing  p  and  $  ).  The  value  of  (  a< t ) |  at  this  point 
is  | p |/4  which  is  at  most  1/16.  Representative  values  of  FQ(t)  for 
xl=1'x2  =1,5  and  p=-25, -.125,0,  .  125,  .25  are  plotted  in  Figure  1. 
The  curve  with p  =0  corresponds  to  the  system  reliability  if  the 
assumption  of  independence  is  true.  Since  most  applications  of 
interest  involve  reliabilities  of  .75  or  greater,  we  plot  in  Figure  2 
the  ratio  of  the  100  pth  upper  percentiles  under  dependence  and 
independence  vs.  the  correlation.  From  Figure  2,  it  appears  that. 
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when  the  predicted  system  reliability  under  independence  is 
greater  than  . 90/  misspecifying  the  dependence  parameter  has 
little  effect.  However,  in  the  range  where  the  predicted  system 
reliability  under  independence  is  less  than  . 75,  misspecifying 
the  dependence  parameter  may  lead  to  errors  exceeding  6%. 

Maximum  values  of  { A  ( t )  |  are  tabled  in  Table  1,  for  X^=l  and 
various  values  of  X2> 

The  mean  time  to  system  failure  based  on  (2.1),  assuming 
independence,  is 

Uj  =  1/(X1+X2)  (2.6) 

and  that  based  on  (2.4)  is 

The  amount  of  error  in  modeling  system  mean  life  is 

6p  X^X2 

VD~VI  =  U1+X2)  (2X1+X2)  (X1+2X2) 

_  6pXlX2wI  (2.8) 

(2X1+X2)  (Xx+2X2) 

whose  absolute  value  obviously  increases  as  |  p  )  increases.  If 
xl=x2 '  this  error  reduces  to  which  has  a  maximum  absolute 

value  of  pj/6. 

It  is  apparent  from  Table  1  and  equations  (2.5)  and  (2.8) 
that  the  error  in  modeling  system  reliability  and  mean  system 
life,  based  on  independence,  increases  as  |p|  increases  and  is  a 
function  of  the  relative  sizes  of  ^  an<*  In  particular,  when 

the  mean  life  of  one  component  is  substantially  greater  than  the 


(Xi+x2) 


^(XJ.H 


(2X1+X2 


(X1+2X2) 


(2.7) 


I 


I 

i 
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mean  life  of  the  second  component,  then  the  behavior  of  the 
system  is  well-approximated  by  the  behavior  of  the  shorter- 
lived  component  acting  alone.  This  can  be  seen  in  (2.4)  and 
(2.7)  by  letting  0.  In  this  instance,  we  see  also,  from  (2.5) 

and  (2.8),  that  the  amount  of  error  incurred  by  assuming  indepen¬ 
dence  is  negligible. 

3.  ESTIMATING  SYSTEM  RELIABILITY  FROM  COMPONENT  DATA 


A  common  practice  in  predicting  system  reliability  is  to 
test  each  of  the  components  independently  and  then  to  use  the 
data  to  obtain  confidence  bounds  on  system  reliability.  These 
bounds,  obtained  by  Mann  and  Grubbs  (1974),  assume  that  the 
component  lifetimes  are  exponential,  and  that  the  components 
act  independently  when  linked  in  series.  In  the  bivariate  case 
the  bounds  are  computed  as  follows:  for  the  j  component  suppose 
that  prototypes  have  been  tested  until  (r^n^)  failures  occur. 


Let  Zj  be  the  total  time  on  test  for  the  jth  component.  Define 
M*=£(rj-1)/Zj+  ^(r.-l/Z^/^rj-D/Z^],  (3.1) 

V*=£(rj-1)/Z?+  £<£(r.-lYzl'}/{£(r.-l)/Z^.  (3.2) 

An  approximate  r  level  lower  confidence  level  for  system  reliability 

at  time  t„  is 
m 

exp[-tm  M*  £l-vV<9M*2)  +  nY(V*)V<3M*)^  3 Jr  (3.3) 

where  n  is  the  100  Y  percentile  of  a  standard  normal  random  variable. 


Table  1.  Maximum  values  of  |  a  ( t )  |  for  x^=l  and  various 
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When  the  system  being  evaluated  has  dependent  components  those 
bounds  may  be  misleading.  The  problem  is  that  component  data  is 
independent,  since  the  components  are  tested  separately,  but  when 
put  together  into  a  system  some  interdependence  may  develop,  of  course 
such  dependence  is  not  detectable,  in  the  absence  of  some  system 
data,  since  the  data  on  components  we  see  are  independent.  To  study 
the  performance  of  the  bound  (3.3)  when  the  correct  system  model  is 
the  Gumbel  model  (2.2)  a  simulation  study  was  performed.  For  each 
simulated  sample  observations  from  exponential  populations  with 
mean  l/\^,  j=i,2,  were  simulated.  The  two  samples  were  generated  inde¬ 
pendently.  The  confidence  bound  (3.3)  was  obtained.  This  was  then 
compared  to  the  true  system  reliability  at  various  p's  obtained 
from  (2.4).  Ten  thousand  such  bounds  were  simulated  for  each  set 
of  parameter  values.  The  results  for  n1=n2=3, 5 , 10 , x1=l . 0 , x2=l . 5 , 
at  tm=0.1  are  reported  in  Table  2.  Here  the  system  reliability 
ranges  from  .7684  at  p  =  -.25  to  .7891  at  p=.25  with  a  value  of  .7788 
at  P=0 . 

The  results  in  Table  2  show  that  at  high  negative  correlations 
the  coverage  probabilities  are  significantly  lower  than  claimed 
under  independence,  while  for  a  high  positive  correlation  the  intervals 

are  conservative.  This  trend  becomes  more  exaggerated  as  n^,n2 
increase.  This  is  due  to  the  fact  that  as  n^ , n2  increases  the 
bound  approaches  the  reliability  under  independence.  As  seen 
in  Section  1,  the  true  reliability  at  t  is  an  increasing 
function  of  p  so  that  asymptotically  coverage  probabilities 
approach  0(orl)  for  p<0(>0).  On  the  practical  side,  while  the 


Jt 
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estimated  coverage  probabilities  for  are  statistically 

significantly  lower  than  expected  for  small  sample  sizes,  they 
are  not  of  sufficient  magnitude  to  cause  great  concern,  especially 
at  Y  =. 95. 

4.  ESTIMATING  COMPONENT  PARAMETERS 

In  this  section  we  are  interested  in  examining  how  the 
independence  assumption  affects  the  magnitude  of  the  estimation 
error  in  estimating  component  reliability  and  mean  life  from 
data  collected  on  series  systems.  That  is,  for  each  system 
put  on  test  we  observe  its  failure  time  and  an  indicator  variable 
which  tells  us  which  component  caused  the  system  to  fail.  We 
are  interested  in  how  varying  degrees  of  dependence  affect  the 
bias  and  mean  squared  error  of  either  parametric  or  nonparametric 
estimators  of  component  survival  and  mean  life  which  were  obtained 
by  assuming  independent  component  lifetimes. 

4 . 1  Parametric  Approach 

Consider,  first,  the  problem  of  parameter  estimation  in  a 
two  component  series  system.  We  assume  the  two  components' 
survival  functions  are  f\(t)  =  exp(-X^t) ,  i  =  1,2,  and  a  life 
test  is  conducted  by  putting  n  systems  on  test.  We  observe  n^ 
systems  to  fail  due  to  failure  of  the  ith,  i  =  1,2,  component. 

Let  T  denote  the  sum  of  all  n  failure  times.  From  Moeschberger 
and  David  (1971)  the  maximum  likelihood  estimator  of  x^  is, 
assuming  independence, 

X ^  =  r^/T,  i  =  1,2 

so  the  estimator  of  component  mean  life,  p^«x^,  is 
ui  =  H/n^  if  ni>0.  (4  .1) 


Now  suppose  that  we  are  in  fact  sampling  from  the  Gumbel 
distribution  (2.3).  For  this  model,  component  mean  life  is 
the  same  as  in  the  independent  case.  The  random  variables  (n^.T) 
are  independent  (the  conditional  distribution  of  T  given  n^  is 
free  of  n^)  and  is  binominal  with  parameters  n  and 
p^  =  P(min(X^,X2)  =  X^) .  For  this  model 


model 

4P  a1-x2)  x  2 


Pi  =  P(X1<X2)  X1{X1+X2  +  (X1+X2) (2XX+X2) (Xx+2X2) }  (4.2) 

with  p2=l-p^.  From  Mendenhall  and  Lehman  (1960)  approximations 
to  the  moments  of  1/n^,  conditional  on  n^>0,  are 


Ed/nJn^O,  -  KTS^IT' 


E  ( 1/n.2  |  n  ■  >0)  =  Jn.~2)  {.n~3)- 
n  (a-1) (a-2) 


(4.3) 


(4.4) 


where  a  =  (n>l)p., 

The  expected  value  of  T  is  given  by  npD  where  pd  is  given  by 


(2.7)  and 


E(T2)  = 


...  „  f~  2+10  p 

-  n[Jv 


8p((2X1+X2)2  +  (X 


1+Tx2)2j]+  n(n“1)WD 


.  (4.5) 


Thus  the  bias  and  MSE  of  p^,  conditional  on  n^>0,  under  this  model 


B(u,)  =  E ( u . -p . )  = 


(n-2) p. 


i  Mi'  Qn-l)p.-l)J  wi' 
-  2g.  (n-2)  pn 


and  MSE(ui)  =  E(T^)  E(l/nf|ni>0)  -  *v± 


(4.6) 


(4.7) 


for  i  =  1,2. 
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We  note  that  for  large  samples 


lim  B(ui) 
n  -*  “  1 


(4  .8) 


lim  MSE(ui)  =  (  lim  B^))2  (4.9) 

n  -*■  <*>  n  -*■  oo 

for  i  =  1,2. 

For  from  (4.6),  we  see  that 


B(y.  )« 


_  1+2 (n-2 ) p/3 


Cn-3) 


y  x  2  ( n  -  2 )  p  y  ^ 

n^3  +  3 (n-3 ) 


(4.10) 


A  similar  expression  holds  for  B  ( u  2 )  •  Nate  that  (4.10)  consists 
of  two  terms.  The  first  term,  reflecting  sampling  error,  is 
positive  for  all  n  and  dominates  the  bias  expression  for  small 
n.  The  second  term,  reflecting  modeling  error,  takes  on  the 
same  sign  as  the  correlation  and  dominates  for  large  n,  approaching 
the  limit  of  2py^/3. 

For  5*  \ 2  >  the  bias  can  again  be  expressed  as  the  sum  of  a 
term  reflecting  sampling  error  and  a  term  reflecting  modeling  error. 
The  sampling  error  depends  upon  the  expected  number  of  component  1 
failures,  np1.  From  (4.2)  we  note  that  p^.  increases  as  increases 
if  X^>X2  and  decreases  if  ^2.<^2*  T^e  moc^el^n9  error  is  reflected  in 
the  asymptotic  bias  which  is  an  increasing  function  of  p.  A  typical 
representation  of  the  bias  of  y^,  as  a  function  of  sample  size  and 
p,  is  in  Figure  3  where  X^=l,  \2=1'5* 


Comments  similar  to  those  pertaining  to  the  bias  may  be  made 
for  the  mean  squared  error  of  y^.  Figure  4  depicts  the  ratio  of 
the  MSE  (ujjp)  and  MSE  (y^|p=0)  as  a  function  of  sample  size  and 


for  x-1, >2-1.5. 
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4  .2  Nonparametric  Approach 

A  second  approach  to  the  problem  of  estimating  component 

parameters  is  via  the  nonparametric  estimator  proposed  by 
Kaplan  and  Meier  (1958).  The  product  limit  estimator,  assuming 
independent  risks,  is  constructed  as  follows.  As  before 
suppose  n  systems  are  put  on  test  at  time  0  and  ni  systems 
fail  due  to  failure  of  component  i.  Let  Xi ( ± ) , • * • , (n^) 


denote  the  ordered  times  at  which  these  n^  events  occur 

and  let  r.,,  . ..,  r.  be  the  ranks  of  those  ordered  survival 
ll 

times  among  all  n  ordered  lifetimes.  The  component  relia¬ 
bility  for  the  ith  component  at  time  x  may  now  be  estimated 
by  the  product  of  the  individual  conditional  survival 
probabilities,  namely  by 

Fi<x>  *  f1  i£  *<xi(l) 


<  j (i,x) 

where  j  ( i ,  x)  is  the  largest  value  of  j  for  which  x^ ^  <x. 

A  special  note  is  needed  to  cover  the  case  of  x^n  ^  not 

being  the  largest  observed  death  or  removal  time.  To  avoid 
this  problem  we  shall  define  F\  (x)  =  0  for  x  greater  than 
the  largest  observed  failure  time. 

If  the  component  lifetimes,  in  fact,  follow  the  Gumbel 
bivariate  exponential,  Klein  and  Moeschberger  (1983)  have 
shown  that  the  Kaplan-Meier  estimator  is  not  consistent. 

For  i  =  1,  the  Kaplan-Meier  estimator  is  not  estimating 
F^(t)  but  rather  another  survival  function,  H^(t)  given  by 


'X2U 


-Xiu 


t  L^p(l-e  *  HWylil  d  u) ,  t>0  (4.11) 


HL(t)  exp{-Ax/  [j.+4p  (l-e"Xlu)  (l-e~X2ufJ 


Note  that  equation  (4.11)  simplifies,  if  X^x^v,  to 


HL(t)  =  e~'pt  [l+4p  (l-e~tt:)2]h  (4.12) 

% 

which  is  increasing  in  o  .  Similarly,  F^ftjis  actually 
estimating  ^(t),  which  is  defined  analogously. 

% 

Measures  of  the  error  in  estimating  f\(t)  by  F^(t)  are 

again  the  bias  and  mean  squared  error  of  §v  (t)  computed 
under  the  dependence  model.  Under  this  model  the  Kaplan - 
Meier  estimator  is  equivalent  to  the  estimator  one  would 
obtain  based  on  n  observations  from  an  independent  system 
with  component  survival  distributions  given  by  (4 .11) 
or,  if  X^=X_2,  by  (4.12).  Hence  from  Kaplan  and  Meier  (1958) 

'X, 

it  follows  that  the  variance  of  F^t)  is  given  by 


5,  t|dH.(u)j 

V(Fi(t))  .  Bpt,2  '  „  gi<u)2 


(4.13) 


Thus  the  bias  and  MSE  of  F^t)  are,  from  (4.11)  and  4  .13) 


B(Fi(t) )=Hi(t)  -  Fi(t) ,  t  >  0, 


(4.14; 


and  MSE(F.  (t) )  =  (H,  (t)  -F.  (t)  )  2  +  H,  (t)  2  tldHi  (u)l  ,t>0. 

i  x  j  — - - 5—  (4. 


0  n  Hj.  (u) 


15) 


The  estimator  is  not  consistent  since  B(Fi(t))  is 
independent  of  n,  and  not  necessarily  zero.  Also 

"X, 

MSE ( F^ ( t ) )  consists  of  a  factor  which  depends  only  on  the  model 
error  and  is  free  of  sample  size,  and,  a  term  which  tends  to  0 
as  n  tends  to  infinity. 


o<  +X0 


Note  that  in  the  case  of  equal  component  lifetimes, 
\^=A2=<*1  the  bias  determined  from  (4.12)  and  (4  .14)  simplifies 


In  the  general  case,  the  integral  in  (3.11)  needs  to  be 
evaluated  numerically.  The  bias  of  the  Kaplan-Meier  estimator 
was  calculated  for  various  values  of  A^  and  p.  A  representative 
plot  of  the  bias  appears  in  Figure  5,  where  \^=1 , Aj^l • 5 ,  and 
| p { =0 ,  .125,  .250.  It  is  apparent  that  the  bias  is  largest 
for  values  of  t  in  the  neighborhood  of  an  interval  which  captures 
the  mean  component  lifetimes.  The  absolute  value  of  the  bias 
ranges  from  0  to  .11,  in  this  example. 

MSE(E\(t))  was  calculated  for  various  values  of  A^,  n,  and  p. 
Its  magnitude  is  typified  in  Figures  6a  and  6b,  where  A1=1,A2=1.5 
and  n=10,  50,  respectively.  For  A^=1,A2=1.5,  and  n=<»  MSE(F  (t)  ) 
may  be  obtained  by  squaring  (4.14)  or  by  squaring  the  ordinate 
values  in  Figure  5.  The  mean  squared  error  of  the  Kaplan-Meier 
estimator  may  be  quite  large  for  small  sample  size  n  and  moderately 
large  for  "large"  p,  the  former  being  a  more  crucial  factor  than 
the  later. 
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ABSTRACT 


A  common  assumption  made  in  modeling  and  analyzing  data  from 
series  systems  is  that  the  component  lives  are  statistically 
independent.  This  study  investigates  the  magnitude  of  the 
errors  one  may  incur  by  erroneously  assuming  the  component  life¬ 
times  have  independent  exponential  distributions  when  in  fact 
the  lifetimes  follow  the  bivariate  exponential  distribution  of 
Marshall  and  Olkin  (1967). 

Key  Words:  System  Reliability;  Competing  Risks;  Marshall- 
Olkin  Bivariate  Exponential;  Series  System. 


1. 


INTRODUCTION 


Consider  a  two  component  system  where  the  components  are 
linked  in  series.  Such  a  system  will  function 
properly  provided  both  components  function  properly.  Suppose 
that,  if  the  two  components  are  tested  separately,  the 
respective  times  to  failure  are,  X^ ,  X£ .  A  common  assumption 
made  in  life  testing  is  that  these  component  failure  times 
follow  exponential  distributions  (cf.  Epstein  and  Sobel  (1954), 
Boardman  and  Kendall  (1970),  Mann  (1974),  Chao  (1981),  etc.). 

Often,  when  modeling  a  system  linking  two  components  in 
series,  the  assumption  that  X^,  X£  are  statistically  independent 
is  invoked.  In  general  such  an  assumption  is  not  testable  due 
to  the  identif iability  dilemma  (see  Basu  and  Klein  (1983)  for 
a  discussion  and  references).  This  assumption  is  important  in 
both  theoretical  modeling  of  the  resulting  system  and  in 
obtaining  estimators  of  component  parameters  from  data  collected 
on  series  systems. 

The  purpose  of  this  study  is  to  investigate  the  effects 
of  departures  from  this  independence  assumption  on  both  the 
modeling  and  estimation  problem.  The  specific  form  of 
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departure  from  independence  which  we  shall  assume  is  that  the 
joint  distribution  of  X^,  X2  is  the  bivariate  exponential 
distribution  of  Marshall  and  Olkin  (1967).  This  distribution  has 
joint  survival  function 

P(Xi>xi,  X2>x2)  =  exp(-.\^x1-X2x2~^, ,  maxtx^.x,,)) 

for  Xv  \2  >  0,  A12  >  0,  xx,  x2  >  0.  (1.1) 

The  marginal  distributions  of  and  X2  are  exponential  with 
parameters  3^  =  +  X^2,  i  =  1,  2.  This  distribution  can  be 

derived  from  both  a  fatal  and  non-fatal  shock  model.  It  has  a 
singular  component  when  is  nonzero  which  reflects  the 
possibility  of  the  system  receiving  a  shock  of  sufficient 
magnitude  to  simultaneously  destroy  both  components.  The 
distribution  reduces  to  that  of  independent  exponentials  when 

A12  =  °* 

In  Section  2  we  examine  the  effects  of  assuming  the 
components  to  have  independent  life-lengths  in  modeling  series 
systems,  on  tne  system  reliability,  system  mean  life  and  crude 
system  life,  when  in  fact  the  distribution  of  lifetimes  follows 

the  bivariate  exponential  distribution  (1.1).  In  Section  3  we 
investigate  the  effects  of  the  independence  assumption  on 
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estimating  component  parameters  when  in  fact  the  bivariate 
exponential  is  the  correct  model. 

2.  MODELING  SERIES  SYSTEMS 

2.1  System  Reliability 

Suppose  that,  based  on  testing  each  component  separately, 
an  investigator  knows  that  the  marginal  distributions  of  the 
component  lifetimes  are  given  by  P(X^  >  x)  =  exp  (-(A.+Al2)x), 
1*1,2.'  Also,  suppose  the  investigator  believes,  for  the 
purpose  of  modeling  system  life,  that  the  components  are 
statistically  independent  when  in  fact  the  joint  distribution 
is  the  bivariate  exponential  (1.1).  Under  the  assumption  of 
independence  the  system  reliability  at  a  mission  time  t  is 
give  by 

P(min(X^,  X2)  >  t  |  Independence)  =  Fj(t)  =  exp(-('*+X^2)t) , 

where  A  =  while  under  the  bivariate  exponential  model 

we  have 

P(min(X^,  X2)  >  t  |  BVE)  -  FQ(t)  *  exp(-At). 

The  amount  of  error  made  in  assuming  independence  is 

FQ(t)  -  Fj(t)  -  exp(-*t)  (l-exp(-.\12t))  (2.1.1) 


which  is  always  non-negacive  so  that  the  system  is  more 
reliable  then  one  would  expect  under  independence. 

For  A^,  A2,  A^  fixed  the  maximum  difference  between 
Fp(t)  and  F^t)  occurs  at  t  =  A“*  lnU+A^/A]  and  the 
magnitude  of  this  difference  is  C  (l+c)-^p/ (1+p)  where 
0  =  ^12^  is  the  correlation  between  X^,  and  X2 .  This  maximum 
difference  approaches  1/4  as  P  increases  to  one. 

Alternatively,  an  investigator  modeling  system  life, 
based  on  component  knowledge,  may  wish  to  assess  the  error 
(2.1.1)  at  some  fixed  mission  time.  Consequently,  for  \^,  >>2, 
and  t  >  0  fixed,  the  maximum  deviation  between  FD(t)  and  f^Ct) 
occurs  at  A^  *(ln2)/t  and  this  maximum  difference  is 
.25  exp(-(A1+A2)t) .  In  the  range  of  values  of  t  where  Fj.(t) 
is  less  than  0.5  the  magnitude  of  this  difference  is  F^(t) 
itself  so  that  the  true  system  reliability  is  twice  as  large 
as  that  predicted  u'ing  the  erroneous  assumption  of  independence. 

2.2  System  Mean  Life 

Under  the  assumption  of  independent  component  lifetimes 
the  mean  time  to  failure  of  the  series  system  is  y  =  (A+A^)"1 
while  under  the  correct  bivariate  exponential  model  the  mean 
system  life  is  yD“  A-1 .  Clearly  u  -y  =  A^/MA+A^)  =  PMj 


which  is  always  non-negative.  The  largest  difference  between 
Up  and  Uj  occurs  when  A^2  =  ( X 2 )  ^ an<^  c^e  ma8niCu^e  of  this 
difference  is  (3-2*2) / (A^+A^)  =  - 1 71 6/( A ^h-A ^ ) .  Again  as  in  the 
case  of  system  reliability  one  may  be  appreciably  misled  when 
P  is  not  zero. 

2.3  P(XX  <  X2) 

Investigators  working  with  series  systems  would  like  to 
know  the  value  of  P(X^  <  X2)  so  that  they  can  determine  the 
dominant  cause  of  system  failure  and  hence  concentrate  their 
efforts  on  improving  the  reliability  of  that  component.  Under 
the  assumption  of  independence  we  have 

*  P(X^  <  X2 1  Independence)  =  (A^+A,^)  /  (A+A^)  while  under  the 
dependence  model  we  have  =  P(X1  <  X2|BVE)  =  A^/A.  The 
difference  between  Trj  and  ^  is  A^2 (A-A^^) /(A(A+A12) )  = 

/(1+P)  which  is  an  increasing  non-negative  function  of 
p  with  a  supremum  of  1/2  as  p  goes  to  unity.  Similar  results 
hold  for  PCXp  >  X2). 

2.4  Crude  System  Life 

An  observable  quantity  in  series  systems  is  the  crude 
system  life  given  by  (t)  =  POninCX^  X2)  >  t,  min(X2 ,X2)=Xj ) , 
j  *  1,  2.  For  the  jth  component  this  is  the  observable 
survival  function  for  systems  where  the  jth  component  fails 
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first.  Under  the  assumption  of  independence  we  have 
Qj  (t !  Independence)  =  Q^(t)  =  ((A^+A^)  /  (A+A^)  )exp(-(A+A12)  t) 
while  under  the  correct  bivariate  exponential  model  we  have 
Qj(t|BVE)  =  QjD^(t)  =  (X^/A)exp(-^t) .  Hence  we  have 


QjI\t)-Q^D)(t)  =  exp  (-At)  ^  exp(-A12t)}. 


For  A^^  not  zercv  the  crude  survival  curve  under  independence 
crosses  the  curve  computed  under  dependence  once  from  above  at 
Cq  *  -A^2  In f A^ (A+A^2) /(A(A^+A12)) ] .  In  the  range  t  <  tQ  the 
maximum  difference  between  1 ^ ( t )  and  QjD^(t)  occurs  at  the 

origin  and  is  the  same  as  the  difference  between  the  respective 
probabilities  that  the  jth  component  fails  first.  In  the 
range  t  >  t^  the  maximum  deviation  between  the  two  curves  is 
at  t^  =  ~A^2  ln[Aj/(Aj+A^2) ]  and  this  maximum  deviation  is 


Xj  X/A12  A12  Xj  m  £_Vj. 

Aj+A^  A  ( A+A^j)  (i+p) 


/ _ 2i_)i/p 

S)j+P 


where  7T  is  the  probability  the  jth  component  fails  first  under 
Che  bivariate  exponential  model. 


3. 


ESTIMATING  COMPONENT  RELIABILITY 


A  second  area  where  Che  assumption  of  independence  plays  a 
key  role  is  in  estimating  component  parameters  from  data  collected 
on  series  systems.  Continuing  in  the  same  vein  as  in  Section  2, 
we  assume  the  components  have  exponentially  distributed  lifetimes 
with  unknown  failure  rates  1  ^ •  We  put  n  systems  on  test 

and  for  each  system  we  observe  the  failure  time,  t^,  and  an 
indicator  variable  which  tells  us  which  component  caused  the 
system  failures. 

If  the  data  is  from  systems  whose  components  are  independent 
then  the  sufficient  statistics  for  are  T  =  It^  and  the 
number  of  systems  where  the  first  component  failed  first.  The 

/s 

maximum  likelihood  estimator  of  0^  is  0^  =  M^/T  and 

3>2  ~  (n-M^)/T.  David  and  Moeschberger  (1978)  show  that  the 

'X)  S' 

estimators  0.  *  ((n-l)/n)0j  are  the  best  unbiased  estimators  of 

% 

0j  and  that  for  n  >  2  the  variance  of  0^  is 


Now  suppose  an  investigator  carrys  out  the  above  analysis, 
believing  the  component  lifetimes  are  independent,  when  in  fact 
the  joint  distribution  of  the  component  lifetimes  is  bivariate 
exponential.  If  the  data  is  from  a  bivariate  exponential 
series  system  there  is  a  positive  probability  that  the  two 
components  failed  simultaneously.  In  such  a  case,  suppose 
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the  investigator,  due  to  his  strong  belief  in  the  independence 
assumption  where  such  simultaneous  failures  are  impossible, 
records  this  occurrence  as  a  failure  from  component  1  with 


probability  $  and  a  failure  from  component  2  with  a  probability 
1— <J> .  Such  action  is  done  independently  of  the  system  under 
consideration  and  this  probability  is  constant  from  one  system 
to  another.  In  this  setting  the  parameters  of  interest  are 
0.  *  ^j+^2’  j  ■  1,  2.  For  the  ith  system,  let  X^,  X2i  be  the 
unobservable  component  lifetimes.  The  recorded  statistics  are 
now  T  =  I  min(X^,  X2^)  and  X2^)  where 


^(Xli>  X2i^ 


1  if  Xu  <  X2i  or  Xu  *  X2i  and  Y.  *  1 
0  otherwise 


where  the  Y^'s  are  independent  and  identically  distributed 

Bernoulli  random  variable  with  PCY^  ®  1)  =  4>. 

Under  the  bivariate  exponential  distribution  T  has  a  gamma 

distribution  with  parameters  n  and  X  so  E(T  ')  *  X/(n-l) , 

E(T-2)  -  X2/ [ (n-1)  (n-2)  ]  if  n  >  2.  Also  EQ^)  «  n (X^X.^) /X 

and  E(M*)  ■  nCX^+^X^) /X+(n-l)  n [X^+^X^J 2/X2 •  Also  T  is 

% 

independent  of  Thus  the  bias  of  0^ ,  under  this  model  is 

B(^)  -  E^j-Gp  *  X12(l -<P)  and  B(&2)  “  ^i2^"  *n  ^oth  cas®s 
the  bias  is  an  increasing  function  of  the  correlation  p  so 
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that  the  more  correlated  the  components  the  greater  the  magnitude 
of  the  estimation  error.  The  respective  mean  squared  errors  are 

MSE(0j_)  *  [(n-l)A(A1+4>X12)  +  (A1+4)A12)2]/[n(n-2)]+(l-<!))2Aj2 

and 

HSE(92)  -  t(n-l)A(A2+(l-0)A12)  +  (A2+(i-0)X12)2]/[n(n-2)]-(42A22 

which  are  quadratic  functions  of  p.  As  n  goes  to  infinity  both 
mean  squared  errors  tend  to  nonzero  constants  so  the  estimators 
are  inconsistent. 
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1.  INTRODUCTION 

Competing  risks  arise  in  a  wide  range  of  life  testing 

problems.  Typical  areas  of  application  are  the  study  of 

series  systems  in  the  engineering  sciences  and  biological 

systems  in  the  medical  sciences.  An  important  area  of 

application  is  the  analysis  of  censored  data  where  some 

systems  or  individuals  are  lost  or  withdrawn  from  a  study 

prior  to  observing  the  endpoint  of  interest.  Competing 

risks  are  often  modeled  by  a  vector  T  =  (T, ,  ....  T„)  of 

-v  x  ? 

nonnegative  random  variables  representing  the-  potential 
times  to  failure  from  each  of  the  p  causes.  We  cannot 
observe  T  directly  but  instead  we  see  the  system  failure 
tine  Y  =  min(T^,  i  =  1,  ...,  p)  and  the  failure  pattern 
£(£)  =  I  such  that  Y  =  for  iel  and  Y  <  if I,  where 
lex  the  set  of  all  subsets  of  1,  ...,  p  .  Based  cr.  this 
information  we  wish  to  estimate  the  marginal  survival 
probabilities  Sj(t)  =  PCminCT,.  ,  iel)  >  t),  t  >  0,  lex. 

A  common  assumption  made  in  analyzing  competing  risk 
experiments  is  that  the  T^'s  are  independent  random 
variables.  Such  an  assumption  is  not  testable  due  to  the 
identif iability  dilemma  (see  Basu  (1981)).  Under  an 
assumption  of  independent  risks  a  consistent  estimator 
of  Sj(t)  is  the  Kaplan  and  Meier  (1953)  product  limit 
estimator.  In  Section  2  we  show  that,  if  the  risks  are 
dependent,  then  the  product  limit  estimator  may  be  an 
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inconsistent  estimator  of  S,(t).  The  quantity  to  which 
this  estimator  converges  is  obtained  sc  that  one  may 
investigate  the  estimator's  robustness  to  departures 
from  independence.  In  Section  3  we  illustrate  such 
robustness  considerations  for  some  well-known  bivariate 
exponential  distributions. 

2.  INCONSISTENCY  OF  THE  PRODUCT  LIMIT  ESTIMATOR 

The  Kaplan  and  Meier  product  limit  estimator  is 
constructed  as  follows.  Let  0  =  Y^^  ±  —  •••  —  ^(n) 

denote  the  ordered  system  failure  times  of  n  systems  put  on 
test.  The  product-limit  estimator  of  STCr)  is 

SjCt)  =  niC(n-i)/(n-i+l)]  (2.1) 

where  the  product  is  over  the  ranks  i  of  those  ordered 
observations  1  <  i  <  n,  such  that  YJ(i)  <  t  and  Y(i) 

corresponds  to  a  death  from  the  simultaneous  cause(s)  j  s  J, 

a 

J  n  I  i  2.  Sj.(t)  is  undefined  for  t  >  if  the  largest 

failure  time  corresponds  to  causes  in  J  where  J  '  I  =  3 . 

If  the  assumption  of  independence  is  correct  and  the 
crude  probability  functions  defined  by  FCt,  I)  =  ?(Y  >  t, 
5((£)  =  I)  have  no  common  discontinuities  then  Langberg , 
Proschan  and  Quinzi  (1981)  [LPQC198D]  have  shewn  that  the 
product  limit  estimator  is  consistent.  They  also  show  that 
if  the  F(t,I)’s  have  no  common  discontinuities  then  for  a 
very  particular  form  of  dependence  structure  the  product 
limit  estimator  is  consistent.  We  note  in  the  following 
theorem,  that  their  results  can  be  used  to  study  the  robust¬ 
ness  of  the  product  limit  estimator  to  departures  from  in¬ 
dependence  and  that,  in  general,  if  the  risks  are  dependent 
then  the  product  limit  estimator  (2.1)  is  inconsistent. 
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THEOREM  1.  Let  T  =  (T,  ,  ,  . .  ,  T„  be  a  vector  of  r.on- 
negative ‘random  variables  with  system  life  Y  =  min(T^,  T^) 

and  failure  pattern 


SCT1  =  I,  if  Y  =  T. ,  iel  anc  V  <  T- ,  ifl  (2.2) 

2 ,  otherwise 

Define  rCt)  =  P(Y  >  t),  FCt,  I)  =  ?<Y  <  t,CC£l  =  I),  Ter, 
t  >  0,  and  let  a(.F)  =  {t:?Ctl  >  0>  be  the  support  of  F. 

For  Is-T  define  I-  =  {Jel:  J  I  *  0}.  3ased  on  a  random 
sample  of  size  n  let  Sj  n(t)  be  the  product  limit 

estimator  (2.1).  If  the  functions  F(‘>I)  have  r.o  common 
discontinuities  on  [0,  a(?)]  then 


sI,aCt)  - 


J  el. 


3jCt) 


a.  s . 


ca,.3) 


where 

ST(t)  =  H  [F(a)/?(a‘)]exD[-/t(dFc(. ,  0  <  t  <  oCF)  . 

J  a<t  0 

(,2,4) 

where  the  product  is  over  the  set  of  discontinuities  of 
F(*,J)  and  FC(*,J)  is  the  continuous  part  of  FO,J). 

PROOF.  The  proof  follows  directly  by  applying  the 
results  of  Langberg,  Proschan  and  Quinzi  (.1978)  [LPQC1973)] 
and  LPQC1981).  By  Theorem  4.1  of  LPQC1978)  T  =  H  where 
H  is  a  vector  of  (2^-1)  independent  components  indexed 
lexicographically  by  lex  with  P(Hj  >  t)  =  S^Ct)  given  by 
a .  4  > . 
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Let  =  CTli ,  ....  T0,- )  ,  i  =  1,  ...»  n  be  independent  and 

identically  distributed  as  T.  Replace  ?  and  F(*,J)  in  (2,4) 

^  s  „n 

by  their  emoirical  counterparts  ?  (t)  =  >  x { „•  >  t}/n 

1  “  n  “  4  -  1  - 


Fr  (t  ,J)  =  l  <  t,  CCT,.)  =  J}/n,tocbtai.n  3-  „(t).  B 

i*l  L  ”  ''•* 

X (A)  is  the  indicator  function  of  the  set  A. 

By  Theorem  4.7  of  LPQC1981]  5-  Ct)  *  GTCt)  a.s.  for 
J  x  ,n  J 

J  e  x  .  Routine  algebraic  minipulation  shows  that 


ere 


I  »n 


Ct)  = 


n  3 


Js  X . 


J  ,n 


(t )  so  the  result  now  follows 


In  general,  as  seen  in  the  examples  in  the  following 

section,  II  Sr(t)  i  STCt)  sc  that  an  investigator  may  be 
JsXT  J  1 

seriously  misled  by  incorrectly  assuming  that  the  component 
lifetimes  are  independent.  This  has  been  noticed  by  Fisher 
and  Kanarek  (1974)  in  the  problem  of  analyzing  clinical 
trials  with  censored  data.  Theorem  1  allows  an  investigator 
to  quantify  the  effects  of  the  independence  assumption  by 
computing  the  right  hand  side  of  C 2. 3)  for  some  plausible 
dependent  models. 

LPQ( 1981)  have  shown  that  for  a  special  type  of 

A 

dependence  the  estimator  Sj  n(t)  is  consistent,  ’^e  state 
their  result,  without  proof,  as  a  corollary. 


,-X 
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COROLLARY  1.  Assume  that  the  conditions  of  Theorem  1 
% 

hold.  Then  Sj  nCti  *  SjCt)  a,s.  if  and  only  if  the 
following  two  conditions  hold. 

i)  STCa)/SjCa“)  =  F(a)/F(a") ,  a  discontinuity  point  of 

£F(G,J),  where  the  sum.  is 
over  Jef j 
=  1  ,  otherwise. 

ii)  PCminCT^,  iel ’  )>t  | aiin(T^  ,  iel)  =  t) 

=  P(min(T;  ,  iel ' ) >t  j  mir.CT^  ,  isT)>t) 
where  I'  is  the  comoliment  of  I-:  I. 


EXAMPLES 


In  this  section  we  present  some  representation  examples  of  the 
use  of  Theorem  1  in  determining  the  effects  ,  on  estimating  marginal 
survival,  of  the  independence  assumption  for  some  bivariate  expon¬ 
ential  life  distributions.  Let  (T-^Tj)  denote  the  time  to  failure 
from  components  1  and  2,  respectively,  in  a  series  system. 

Let  F(t^,  t2^  be  the  joint  survival  function  of  (T^,  T,) 

and  S^(t)  =  P(T^>t),  i  =  1,  2  the  marginal  survival  functions. 

* 

Let  S^Ct)  be  the  estimator  (2.1)  of  S^(t),  and  let  (t), 

J= {1} , {2 } , (1 ,2 }  be  given  by  (2.4).  Then  S^_(t  ,.  j  (t  )S{  ^ ^a,s 

by  Theorem  1,  if  the  functions  F(t,  I)  have  no  common 
discontinuities.  Note  that  if  P(T^  =  Tj)  =  0  then  5/  i_  ^ t  ^ =  ^ 
for  all  t.  We  now  give  some  examples 


♦ 

*  1 


T - 


EXAMPLE  1.  CSlock  and  3asu  (.1974)) 

Let  ?(t^,  tj)  =  [A/ CA^+A2  )  ]exp(-X1  t^-Aj  t2~A12Inax^  “1 ’^2  ^ 
-  [An/ (A^  ^A^  )  lexpC-AniaxCt,  ,t,  ) )  ,  for 


^  »  ^2  >  ^  ^  1  9  —  ^  ^  =  ^  1  +  ^0  *  A  ,  ,,  . 


Here  S^(t)  = 


by  theorem  1 


12  - 


12 


•1*  w2 


rr-n2)  «?<-oi*»,2>t)  -  rT^|p 


>  C  . 


t  >  0 ,  but 


A,  A 


S^Ct )  -  exp(-Q-  ^r— y  t),  a.  s .  t  >  0,  i  =  1,  2. 

EXAMPLE  2.  G umbel  (1350) 

Let  P(t^,  t2)  =  expC-A^t^  -  A7t2  -  A^t^tj)  t,  ,  t2  >  0, 

A i^* X 2  >  0  ,  A12  >  0  . 

Here  S^t)  =  exp(-Ait),  t  >_  0  but  by  Theorem  1 
2  2 

S^(t)  exp(-  Ait  "  A12  t  /2 ) ,  a.s.  for  t  >  0,  i  =  1,  2. 

EXAMPLE  3.  G umbel  (1960) 

Let  PCt^,  t2)  =  exp(-A1t1-A2t2 )Cl+A12 

-A12 (exp(-A^t^)  +  exp( -X?t2 ) ) ] , 


+  A, 2 exp (-A^t^-Ajtj)] 


tl»  r2  >  °>  H’  *2  *  3’  A12  -  3  ' 


Tzr 


«*  . a  i  i 
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Here  S^Ct)  =  expC-A^t)  but  by  Theorem  1 

-X4t  -A . t 

Cl+A  (l-e  J  )(l-2e  1  )] 

S-Ct)  ■*  exp{-\.  /  - 1 -  dt  where 

1  0  -\±z  -X.t 

r l+A1?(l-e  HI-  e  J  ) 

j  is  the  compliment  of  i  in  {1,  2}. 

EXAMPLE  4,  Marshal-0  IJcin  (1367) 

Let  F(t^,  t2)  =  expC-A^  ti“*21:2  *  X12  max^t,,  1 2 ) )  » 
d  ^  >  tj  *  0>  X  ^ ,  X  ^  0,  ^12  ^  ^  • 

Here  S^Ct)  =  expC-CA^  +  X^Jt).  In  this  case  the  conditions 
of  Corollary  1  are  met  so  S^(t)  ■*  S,-(x)  a.s. 


1 
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by 
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Abstract 

The  problem  of  testing  for  independence  of  the  component  lifetimes 
when  the  components  are  linked  in  series  is  considered.  To  avoid  the 
problem  of  nonidentifiability  the  marginal  component  lifetimes  are  assumed 
to  be  known.  In  this  setting  a  modified  version  of  Kendall's  Tau  is 
proposed.  This  test  statistic  is  obtained  by  replacing  those  component 
lifetimes  which  cannot  be  observed,  due  to  system  failure,  by  conditional 
probabilities  computed  under  independence.  A  small  scale  simulation  study 
of  the  power  of  this  test  shows  the  test  has  reasonable  power  for  relatively 
small  sample  sizes. 


Key  Wbrds:  Series  Systems;  Test  for  Independence;  Kendall's  Tau; 
Ejponential  Distribution. 
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1.  INTRODUCTION 

A  common  assumption  rrade  in  modeling  series  systems  is  that  the 
component  lifetimes  are  statistically  independent.  This  independence 
assumption  is  also  routinely  made  in  analyzing  data  collected  from  series 
systems.  Recently,  Klein  and  Moeschberger  (1983)  and  Moeschberger  and 
Klein  (1982)  have  shown  that  one  may  be  appreciably  misled  by  erroneously 
assuming  the  independence  of  component  lifelengths  in  modeling  and  estimating 
the  parameters  of  certain  bivariate  exponential  series  systems.  Thus  it 
is  desirable  to  have  a  test  of  this  assumption  of  independence  based  on 
data  from  series  systems. 

If  no  assumption  about  the  underlying  distribution  of  the  component 
lifetimes  is  made  it  is  impossible  to  test  for  independence  due  to  the 
identifiability  problem  (see,  e.g.  Tsiatis  (1978),  Miller  (1977)  or 
Basu  and  Klein  (1982)).  That  is,  given  any  set  of  observable  information, 
collected  from  a  series  system  with  dependent  component  lifetimes ,  there 
exists  a  series  system  with  independent  component  lifetimes  which  yields 
the  same  set  of  observable  information.  Hence,  any  test  of  independence 
must  incorporate  some  additional  information  about  the  underlying  model 
to  be  able  to  detect  departures  from  independence. 

In  Section  2  a  modification  of  Kendall's  (1938)  test  for  independence 
is  considered.  This  modification  assumes  that  the  marginal  component  life 
distributions  are  completely  specified.  This  information  could  be  obtained 
by  testing  each  component  separately,  as  if  often  done  in  the  development 
stages  of  system  design  (see,  e.g.  Easterling  and  Prairie  (1971),  Mastran  (1976), 
and  Miyamura  (1982)).  In  Section  3  a  small  scale  simulation  study  shows 
that  the  test  performs  fairly  well  for  relatively  snail  sample  sizes. 


2.  THE  TEST  PROCEDURE 


Suppose  that  n  two  component  series  systems  are  put  on  test.  Let 


X.,  denote  the  potential  (unobservable)  failure  times  of  the  first  and 
second  components  of  the  i^1  system.  We  are  not  allowed  to  observe 
(Xf,  Y^)  directly,  but  instead  we  observe  T^  =  min(X^,  Y^) ,  the  system 

1  if  T^  =  ,  the  cause  of  the  system  failure. 


failure  time  and 


Xi  = 


0  if  T.  =  Yi 


Also  suppose  that  the  marginal  survival  functions  of  X^  and  Y^,  F(x)  =  P(X^>x) 
and  &(y)  =  P(Y^>y),  i  =  l,...,n  are  known. 

If  we  could  observe  both  X^  and  Y^  then  a  test  of  independence,  due 
to  Kendall  (1938),  is  to  count  the  number  of  concordant  pairs  and  the  number 
of  discordant  pairs.  A  pair  (Xp  Yi>,  (Xj,  Yj )  is  concordant  if  Xi  -  X_. 
and  Y^  -  Yj  have  the  same  sign  and  is  discordant  if  these  differences  have 
different  signs.  The  test  statistic  is  then  the  number  of  concordant 
pairs  minus  the  number  of  discordant  pairs. 

If  the  data  comes  from  a  series  system  then  only  T^,  1^  is  observed. 
Suppose  we  consider  a  pair  (T^,  1^),  (Tj ,  I.)  with  T^  <  Tj .  If  1^  =  1  and 
Ij  =  1  then  we  know  that  X^  =  T^  <  Xj  =  Tj ,  and  X^  <  Y^,  Xj  <  Yj .  This 
pair  would  be  concordant ,  regardless  of  the  value  of  Y^ ,  if  T^  <  Y^  <  Tj . 

If  Y^  >  Tj  concordance  or  discordance  depends  on  the  value  of  Yj.  Under  the 
null  hypothesis  of  independence,  the  conditional  probability  that  the 
pair  is  concordant  is  CS(T^)  -  G(Tj )  ]/S(T^) .  When  1^  =  1  and  Ij  =  0  then 
Ti  =  Xi  <  =  Tj»  Xi  <  Y^,  Yj  <  Xj.  Here  if  T^  <  Y^  <  Tj  the  pair  would 

be  concordant  and  if  Y^  >  Tj  the  pair  would  be  discordant ,  whatever  the 
value  of  Xj.  Under  independence  the  conditional  probabilities  of  these 
two  events  are  [GCT^  -  G(Tj)]/G(T^)  and  G(Tj)/G(T^),  respectively.  Should 
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1^=0  similar  probabilities,  involving  F,  could  be  obtained.  This 
motivation  suggests  the  following  score  function  for  Tl  <  T4 


CG(T.)-G(T.)]/G(T.) 

if  I.  =  =  1 

X  j  X 

[F(T. )-F(T. ) ]/F(T. ) 

if  I.  =  I.  =  0 

H 

i_i 

H 

1  3 

[S(Ti)-2^(Tj)]/G(Ti) 

if  I.  =  1,  I.  *  0 

[F(Ti)-2F(Tj)]/F(Ti) 

if  I.  *  0,  I.  *  1 

^i’li’Vj)  = 


and  similarly  for  'll  >  TV . 

The  modified  version  of  Kendall's  test  statistic  is 

I  l  <KTi,Ii,T.,I.)/(5). 

l<i< j  <n  133 


(2.1) 


(2.2) 


In  the  appendix  we  show,  that  under  the  null  hypothesis  of  independence  of 
the  component  lifetimes: 

A)  E(t)  =  0 


F(y)G(y)2dG(y)dF(x) 


G(y)F(y)2dF(y)dG(x) 


r,oo 

G(x)2dF(x)  +  \  F(x)2dG(x) 

00  J  -00 

-  2  I  G(x)'1|F(y)G(y)2dG(y)dF(x) 

J_oo  Jx 

-  2  [  F(x)~1f  G(y)F(y)2dF(y)dG(x) 

■)_»  J  X 

{CD  {CD 

+  4(n-2)  F(x) 3G(x)dG(x)  -  2  F(x)G(x)dG(x) 

j  —00  —00 

{DO  {CD 

-  2  F(x)2G(x)2dG(x)  +  3  G(x)2F(x)dG(x) 

j  —00  _©o 

rl  {CD 

F(x)2dG(x)  -  F(x)3dG(x)}  ,  where  F(x)  =  1  -  F(x) 

_ 00  *  —00 


(2.3) 


-  2  [  F(x)2G(x)2dG(x)  +  3  |  G(x)2F(x)dG(x) 


and  G(x)  =  1  -  G(x). 


The  asymptotic  normality  of  t  follows  by  the  results  of  Hoeffding  (1948). 
Hence,  a  test  of  independence  versus  dependence  rejects  if  j t//kt) |  is 


( 
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greater  than  the  appropriate  percentage  point  of  a  standard  normal  random 
variable.  A  test  of  independence  versus  positive  dependence  rejects  if 
T  /Al  i)  is  too  large. 


The  variance  of  x  (2.3)  can  be  expressed  explicitly  in  several  cases. 
Case  1.  F(x)  =  G(x).  In  this  case  (2.3)  reduces  to 

vc?)  =  •  <2.4) 

Case  2.  (Lehmann  structure)  F(x)  =  G(x)a.  Here  (2.3)  reduces  to 


-i  \,r/b\  _  8a[35ct  +  n(9<x  +2a+9)] 
n(n-l)V(T)  -  3 '( 3a+'i‘5"( 3 4-^ )  ( 2a+  3 )  ( 3a+  2 ) 

. . s  _  -6y  v 


(2.5) 


Case  3.  (X,  Y  exponential),  F(X)  =  e”  ,  G(y)  =  e 


hen  (2.3) 


reduces  to 


n(n-l)V(t) 


8X9C35X0  +  n(9X2+2X0+902) ] 

"  3(3X+0y(X+39) (2X+39) (3X+20)  * 


(2.6) 


3.  SIMULATION  STUDY 


A  simulation  study  was  conducted  on  the  AMDAHL  470  computer  at  The  Ohio 
State  University.  The  study  was  performed  by  generating  1,000  samples  of 
n  =  20  or  40  series  systems  with  exponentially  distributed  component  life  times 
F(x)  =  e  x,  and  G(y)  =  e  0  =  1.,  2.  The  value  of  Kendall's  r  =  0,  .3,  .6, 
was  obtained  by  using  the  procedure  of  Johnson  and  Tenebein  (1981). 

Table  1  summarizes  the  results  of  this  study.  Emperical  powers  times 
1000  of  the  modified  test  of  HQ:  independence  versus  H^:  t  >  0  were 
computed  by  comparing  t///(t)  to  the  upper  percentage  point  of  a  standard 
normal.  The  resulting  emperical  powers  (and  their,  standard  errors  in  (  ) 
along  with  the  estimated  expectation  of  9  (and  its  standard  error  in  (  ))  are 
reported. 


TABLE  1.  SIMULATION  STUDY  OF  T 


e 

n 

T 

a  =  .1 

ct  =  .01 

E(t) 

20 

0.0 

97(9) 

44(7) 

17(4) 

7(3) 

-.0039( .0027) 

1 

40 

0.0 

34(9) 

42(6) 

22(4) 

7(3) 

-.0027( .0034) 

1 

20 

0.3 

356(15) 

241(14) 

168(12) 

86(9) 

. 0794( . 0029) 

I 

40 

0.3 

514(15) 

376(15) 

275(14) 

177(12) 

.0766( .0021) 

1 

20 

0.6 

726(14) 

593(16) 

461(16) 

317(15) 

.1596( .0028) 

1 

40 

0.6 

902(9) 

831(12) 

734(14) 

591(15) 

.1596( .0020) 

1 

20 

0.9 

880(10) 

822(12) 

733(14) 

612(15) 

•  225K.0029) 

1 

40 

0.9 

998(3) 

970(5) 

951(7) 

899(9) 

. 2214( . 0020) 

2 

20 

0.0 

92(9) 

46(7) 

20(4) 

4(2) 

-.0028( .0027) 

L 

40 

0.0 

101(9) 

55(7) 

23(5) 

9(3) 

. 0009( . 0019) 

2 

20 

0.3 

426(16) 

319(15) 

233(13) 

139(11) 

.0952( .0030) 

2 

40 

0.3 

599(15) 

463(16) 

338(15) 

222(13) 

. 0891( . 0021) 

2 

20 

0.6 

802(13) 

711(14) 

605(16) 

478(16) 

.1923( . 0029) 

n 

i- 

40 

0.6 

964(6) 

929(8) 

875(10) 

784(13) 

•1910( .0020) 

20 

0.9 

998(1) 

997(2) 

986(4) 

950(7) 

. 3154( . 0021) 

-) 

L. 

40 

0.9 

1000(0) 

1000(0) 

1000(0) 

1000(0) 

. 3169( .0014) 

Nets:  Numbers  reported  are  the  estimates  (standard  error  of  the  estimates). 


Several  conclusions  can  be  drawn  from  this  study.  First,  the  normal 

A 

approximation  to  the  null  distribution  of  t  appears  to  hold  quite  well  for 

A 

relatively  small  sample  sizes  like  20.  Secondly,  E(t)  is  positive  for 
t  >  .3  and  in  fact  we  conjecture  this  is  true  for  x  >  0.  However,  as 

/v 

expected  x  is  underestimating  x  by  a  substantial  margin  and  as  such  is  not 
a  good  point  estimator.  Thirdly,  x  has  reasonably  good  power  when  x  is 
large.  lastly,  the  power  at  fixed  n,  and  x  is  higher  when  the  components  are 
not  identical.  This  is  reasonable  since  in  this  case  our  test  is  based 
on  T  =  minimum  of  (X,  Y)  and  not  on  the  distribution  of  (T,  1=1)  and 
(T,  I  =  2)  the  crude  probabilities. 

4.  CONCLUDING  REMARKS 

This  article  presents  a  technique  for  testing  the  independence  of 
component  lifetimes  based  on  data  collected  when  the  components  are  linked 
in  series.  The  problem  of  nonidentif lability  which  prohibits  testing, 
nonparametrically ,  this  hypothesis  of  independence  is  circumvented,  in  most 
cases,  by  assuming  the  marginal  component  life  distributions  are  known. 

In  general,  this  knowledge  of  the  component  distributions  is  not  sufficient 
to  resolve  the  identif lability  dilenna.  That  is,  there  is  a  small  class  of 
distributions,  mathenatically  contrived,  which  have  the  given  set  of  marginal 
distributions  and  lead  to  the  same  observable  information  as  a  set  of 
independent  component  lifetimes.  Our  test  will  have  no  power  to  detect 
these  departures  from  independence.  However,  for  a  given  set  of  marginals 
there  is  a  large  class  of  joint  distributions  where  knowledge  of  the  marginals 
and  the  system  information  is  sufficient  to  determine  independence.  This 
class  contains  the  standard  bivariate  distributions  of  interest  in  modeling 
series  systems.  For  distributions  in  this  class  our  simulation  study  shows 
that  this  test  performs  well  for  relatively  small  sample  sizes. 
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APPENDIX 


To  show  E(t)  =  0  and  V(x)  =  (2.3)  under  Hq  consider  the  pairs  (Tp  1^) , 

(T2>  I2)  and  (T3,  I  ).  Let  <  Tp  ^  =  I2  =  1}, 

Ag  *  {Tx  <  T2,  II  =  1,  I2  =  0>,  Ag  =  {Tx  <  T2,  =  0,  I2  =  0}  and 

A4  =  <  T2,  I,  =  0,  I2  =  1}.  In  terms  of  the  unobservable  component 

lifetimes,  (Xp  Y^),  A^  =  {X^  <  X2>  <  Yp  X2  <  Y2), 

A2  =  {xl  <  Y2>  Xj_  <  Y1,  Y2  <  X2>,  A3  =  {Yx  <  Y2,  Y]_  <  Xp  Y2  <  X2>,  and 

A^  =  {Y^  <  X2,  Y-J^  <  Xp  X2  <  Y2>.  Under  the  null  hypothesis  of  independence 

is  equally  likely  to  be  either  snaller  or  larger  than  T2  so 


\  E(4>(TpIpT2,I2»  =  | 


^(Xj^J-GCxj) 
Aj  SiCXj) 


dF(x^)dF(x2)dG(y1)dG(y2> 


G(x1)-2^(y2) 
G(x^) 

r  r(yi>-?(y2> 
S  Ryp 


1 


F(y1)-2F<x2) 


dE(x1)dF(x2)dG(y1)dG(y2) 
dF(xL)dF(x2)dG(y1)dG(y2) 
dF(x1)dF(x2)dG(y1)dG(y2) . 


(A.l) 


=  J±  +  +  J3  +  J4  (say). 

Now,  consider 

J.  +  J2  =  I  {|DS(x)-5(y)35(y)dF(y)  +  J  [G(x)-2G(y)]F(y)dG(y)}dF(x).  (A.2) 

1  J_»  Jx  Jx 

Integrating  the  first  inter  integral  in  (A.2)  by  parts  yields  the  negative 
of  the  second  inter  integral  so  that  +  J2  =  0.  Similar  computations  show 
that  J3  +  J4  =  0.  Thus  E($(Tp  Ip  T2,  I2))  and  hence  E(x)  are  both  0. 
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Now 

J  E(*2(T1,I1,T2,I2))  = 

1 2 


G(x^)-G(x2> 


*2 


6(x1)-2G(y2) 


G(V 

T(y1)-F(y2) 


nyx) 


'7(y1)-2F(x2) 


\l  F(y2) 


M  G<*1>  J 

dF(x1)dF(x2)dG(y1)dG(y2) 

dF(x1)dF(x2)dG(y1)dG(y2) 

dF(x1)dF(x2)dGCy1)dG(y2) 


dF(x1)dFCx2)dG(y1)dGCy2) 


(A. 3) 


=  J£  +  J'  +  J3  +  J'  (say). 

After  a  little  simplification  we  have 


j;  +  J'  =  f  ±-  {f° C^(x)-G(y)]2G(y)dF(y)  +  [  [G(x)-2G(y)]2F(y)dG(y)}dF(x) 

1  ‘  J-*5(x1)  ^x  ■'x 


=  f  -  tf  G2  (y)F(y)dG(y)  JdFCx) 


(A.4) 


after  integrating  the  first  term  by  parts.  Writing  F(y)  =  1  -  F(y)  it 
follows  that 


G(x)2dF(x) 


.fi  n 


G(y)2F(y)dG(y)dF(x). 


(A.  5) 


’-•^(xVx 

Similarly  is  of  the  form  (A.  5)  with  the  rolls  of  F  and  G  reversed. 

Hence 

EC^CTpIpT^Ij))  =|[  G(x)2dF(x)  +  |  j  F(x)2dG(x) 

*  ^00  *  «00 

«oo  tco  poo 

-  2  f  -  j  F(y)(a(y)2dG(y)dF(x)  -  2  J  i  |G(y)F(y)2dF(y)dG(x) . 
J-«G(x)  'x  J-*F(x)  'x 


(A.  6) 
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+ 


[F(y,)-2F(x1)]  [F(y5)-2F(x, )]  3 

- =2 - - - - - —  n  dF(x.  )dG(y  . ) 

F(y2)  F(y3)  i=l  1  1 

{y2  »y3<x1  »y2<x2  ’*3^3  ’*1^1 > 


f  CF(y0)-2F(x1 )]  [G(x. )-2G(y,)]  3 

+  2  - i - i - ± - n  dF(x.)dG(y.) 

J  F(y2)  G(x1)  i=l  1  1 

{y2<x1<y3 ,y2<x2 ,x1<y2 ,y3<x3 } 


+ 


[F(y,  )-2F(x0)]  CF(y,)-2F(x,)]  3 

- L - £ - L - n  dF(x.)dG(y.) 

F(y1)  F(y1)  i=l  1  1 

^yl<x2 ,x3 ,yl<xl ,x2<y2 ,x3<y3 ^ 


r  [G(x0)-2G(y, )]  [G(x~)-2G(y, )]  3 

+  - - - ± - - - n  dF(x.)dG(y.) 

J  G(x2)  G(x3)  i=l 

{x2 ,x3<y1,y1<x1 ,x2<y2  ,x3<y3) 


r  CG(x9)-2G(y,)]  [F(y, )-2F(x.,)]  3 

+  2  - £; - ± - i - —  n  dF(x.  )dG(y  . ) 

J  G(x2)  F(xx)  i=l  1  1 

{x2<y1<x3 ,x2<y2 ,y1<x1 ,x3<y3 } 


+  2 


[G(x,  )-G(x„)]  [G(x,  )-2S(y~)]  3 

- i - £ - £ - n  dF(x.)dG(y.) 

G(x]_)  G(x]_)  i=l  1  1 


{x1<x2 ,y3 ,x1<y1 ,X2<y2  ,y3<x3} 


+  2 


C^<x2)-5(x1>]  lF(y^-2T(x^)l 

^(Xj)  ?(y3) 

{x2  .y^x^^ ,  x1<y1 ,  x2<y2  ,y  3<x3  > 


3 

n  dF(x.)dG(y.) 
i=l  1  x 


f  [SU,)-^)]  [G(x.  )-2S(y.) ]  3 

♦  2  - - - - - - - =—  n  dFCx.  )dG(y. ) 

J  G(x2)  SCx^  i=l  1  1 

{x2<x1<y3  ,x2<y2  ,x1<y1,y3<x3} 


(A.  14) 


(A. 15) 


(A.  16) 


(A. 17) 


(A. 18) 


(A. 19) 


(A.20) 


(A.  21) 
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[G(x1)-G(x2)]  CF(y3)-2F(x1)] 
G(xl)  F(y3) 

(y3<x1<x2  ,y3<x3,x],<y1,x2<y2  } 

CF(y1)-2F(x2)]  CF(yi)-F(y3)] 
F(yi)  F(Vl) 

{y1<x2>y3,y1<x1,x2<y2,y3<x3} 

f  CG(x2)-2G(y;L)3  [F(y3)-F(y1)] 
J  G(X2>  F(y3) 

{X2  >y  3<Yl >  x2<y 2  >  y  3  <x3 

CG(x2)-2G(Y;l)3  CF(yi)-F(y3)] 
G(x2)  F(y^) 

{x2<yl<y3 ,X2<y2 ,yl<Xl ,y3<x3} 

,  [F(y3)-F(yi)]  [F(y]L)-2F(x2)3 

J  F(y3)  Tr(y1) 

(y3<yl<X2 ,y3<x3 ,yl<xl ’ x2<y2  * 


Now,  by  combining  like  terms  and  a 

that 


3 

II  dF(x.)dG(y. ) 
i=l  1  1 


3 

n  dF(x. )dG(y. ) 
i=l  1 


3 

n  dF(x.  )dS(y • ) 
i=l  1  1 


3 

II  dF(x.  )dG(y. ) 
i=l  1  1 


3 

II  dF(x.  )dG(y. ) . 
i=l  1  1 


integation  by  parts  one  can  verify 


(A. 7)  +  (A.  13)  +  (A.  19)  =  0,  (A. 10)  +  (A.16)  +  (A.23)  =  0, 

(A. 9)  +  (A. 15)  +  (A.  21)  +  (A.  22)  =  0, 

(A.  12)  +  (A, 18)  +  (A.  25)  +  (A. 26)  =  0, 


(A.  8)  +  (A.  14)  +  (A. 20)  =  f  S(x)G2(x)[2F(x)-l32dF(x) 

J  .00 

(A.  11)  +  (A.  17)  +  (A.  24)  =  j  F(x)F2(x)[2G(x)-l32dG(x). 

J  — OO 


(A.22) 

(A.23) 

(A.  24) 

(A.  25) 

(A.  26) 


(A. 27) 
(A. 28) 
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Integrating  (A.  27)  by  parts  and  combining  the  result  with  (A. 28) 
we  have 


E(b(T1,I1,T2,I2)(()(T1,I1,T3,l3))  = 

-CO  r  I* 00  ^ 

£  F(x)3G(x)dG(x)  -  2  F(x)G(x)dG(x)  -  2  F(x)2G(x)2dG(x) 

4  _oo  *  — OO  *  —to 

f°°  7  F  7  F  3 

+  3  G(x;  F(x)dG(x)  +  F(x)ciG(x)  -  F(x)JdG(x). 

—03  J  — OO  J  — OO 


(A. 29) 


To  find  V(t)  =  E(E  Z  <|>(T.  ,1,  ,T. , I. )/(!?)  )2 
i<j  3  3 

=  C2E(<(,(T1,I1,T2,I2)2)  +  4(n-2)E(<))(T1,I1,T2,I2)<j)(T1,I1,T3,l3))]/n(n-l) 
which  is  equal  to  (2.3)  after  making  the  substitutions  (A.6)  and  (A. 29). 


